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If the effective cosmological constant is non-zero, our observable uni- 
■ verse may enter a stage of exponential expansion. In such case, regions of 



0^ 



it may tunnel back to the false vacuum of an inflaton scalar field, and in- 
flation with a high expansion rate may resume in those regions. An "ideal" 



> 
(N 
0^ 
(N 

■ eternal observer would then witness an infinite succession of cycles from 

o 

l> ■ false vacuum to true, and back. Within each cycle, the entire history of a 

• hot universe would be replayed. If there were several minima of the infla- 

Q ■ ton potential, our ideal observer would visit each one of these minima with 

■ a frequency which depends on the shape of the potential. We generalize 
CZ ■ 

the formalism of stochastic inflation to analyze the global structure of the 
5_j ■ universe when this 'recycling' process is taken into account. 



I. INTRODUCTION 

Inflationary models are designed to produce a universe which is sufficiently homoge- 
neous on all observable scales However, on much larger scales the universe is expected 
to be extremely inhomogeneous. The evolution of the field 0, whose vacuum energy drives 
inflation, is influenced by quantum fluctuations. These fluctuations can be pictured as 
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a random walk of (j) superimposed on its slow roll down the slope of its potential. As a 
result, thermalization of the vacuum energy does not occur simultaneously everywhere in 
the universe, and at any time there are parts of the universe that are still inflating 

On very large scales, the universe is expected to consist of isolated thermalized 
regions embedded in the inflating background. The boundaries of the thermalized regions 
expand into this background, and new regions are constantly being formed, but the high 
expansion rate of the intervening inflating domains prevents these regions from filling up 
the universe. Thermalization inevitably occurs at any given co-moving location, and the 
co-moving volume of the inflating regions decreases exponentially with time. At the same 
time, the physical volume of these regions is exponentially growing. The geometry of the 
inflating regions is that of a self-similar fractal of dimension < 3 [Q. It is illustrated in 
Fig.l for the case of "open" inflation, where the false vacuum decay occurs through bubble 
nucleation . For "new" or "chaotic" inflation the picture would be similar, except the 
thermalized regions would have irregular shapes. 

In the present paper we are going to argue that this picture of the superlarge-scale 
structure of the universe can be significantly modified by quantum fluctuations that bring 
localized parts of already thermalized regions, such as our observable universe, back to 
the inflating false-vacuum state. The modification is particularly important in models 
where the post-thermalization true vacuum is characterized by a positive vacuum energy 
(cosmological constant). In this case the thermalized regions asymptotically approach de 
Sitter geometry, and the rate of fluctuations back to the false vacuum (per unit spacetime 
volume) approaches a constant. Even with an exceedingly small rate, the probability 
for true vacuum to survive at any co-moving location is exponentially decreasing with 
time. Hence, (almost) all the co-moving volume of thermalized regions will eventually be 
recycled back to the inflationary phase. Each nucleated false vacuum region will serve as 
a seed for a new eternally inflating domain, whose internal structure will resemble that 
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shown in Fig.l. The thermahzed regions formed in this domain will in turn produce new 
false vacuum seeds, etc. We call this kind of model a recycling universe. 

Quantum nucleation of regions with a higher energy density cannot occur from a flat- 
spacetime vacuum characterized by a vanishing cosmological constant: such processes are 
forbidden by energy conservation. However, upward fluctuations of this kind can occur in 
an expanding cosmological background, and have been previously discussed by a number 
of authors |17|,|^J^,3-[TT| . The most relevant for our purposes here is the paper by Lee and 
Weinberg |Tl|] who considered a model of a scalar field (f) with a potential V{(f)) shown in 
Fig. 2. Note that both false and true vacua have positive energy densities, pf > pt > 0. 
It has been known for some time that the high-energy false vacuum at = can decay 
by nucleation of true vacuum bubbles. The corresponding instanton ( "bounce" ) has been 
found by Coleman and De Luccia |T^. The bubble nucleation rate is given by 



Aexp[-St + Sfl (i; 



where St, is the bounce action and Sf = — 3/8p/ is the action of the Euclideanized false- 
vacuum de Sitter space (we use Planck units throughout the paper). Lee and Weinberg 
conjectured that the same instanton also describes the inverse process of true vacuum 
decay, where false vacuum bubbles nucleate in a true vacuum background. The nucleation 
rate suggested by Eq.(|lD is 



Aexp[-St + St] (2) 



with St = — 3/8pi. Lee and Weinberg argued that the pre-exponential factors in Eqs. ([1|) 
and (0) are the same. These conjectures were later verified [0 in the case of (1 + 1)- 
dimensional universes, where bubble nucleation can be identified with the production of 
particle-antiparticle pairs. Note that the rate (H) vanishes if the true vacuum has zero 
energy. 
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As it stands, the potential in Fig. 2 is not suitable for inflationary cosmology. This 
potential has no slow-roll region, so most of the vacuum energy remains in domain walls 
and never gets thermalized. We shall see, however, that models incorporating both re- 
alistic inflation and true vacuum decay can be constructed by a trivial modification of 
"open" inflationary models. Moreover, we shall argue that nucleation of inflating regions 
is possible even with the simplest slow-roll potentials, for which the Coleman-de Luccia 
instanton does not exist. 

The recycling nature of inflationary universe may have important implications for the 
question of whether or not the universe had a beginning in time. As we already mentioned, 
inflation is generically eternal to the future, so it is natural to ask if the inflationary 
models can be continued into the infinite past, resulting in a "steady-state" non-singular 
cosmology. This possibility was discussed in the early 80 's, soon after the inflationary 
scenario was proposed, with the conclusion that the idea could not be implemented in the 
simplest model in which the inflating universe is described by an exact de Sitter space 
IT^J^. A more general proof of impossibility of steady state inflation was given in Refs. 
p!5|-[T7|, but we shall see later in this paper that some of the assumptions made in the 
proof do not apply in the case of a recycling universe. The question of the necessity of 
the beginning is therefore re-opened. 

Recycling may also be relevant to the question of making predictions in an inflation- 
ary universe. Recently, there have been a number of attempts [p!8|-p3[] to find probability 
distributions for cosmological parameters such as the effective cosmological constant A or 
the density parameter Q. These "predictions" are based on the principle of mediocrity 
|]20|j2^ , by which we are most likely to live in the most abundant type of civilization that 
can result from the thermalization of a fase vacuum. However, in the inflationary uni- 
verse, there will be an infinite number of infinite thermalized regions, and one faces the 



difficulty of comparing infinities |ll9| , p!8| , p7| . Regularization procedures were introduced 
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in to deal with this problem, but these cannot be directly applied to a recycling 

universe. 

The paper is organized as follows. In the next Section we shall give some examples 
of inflationary models which allow nucleation of false vacuum bubbles. The geometry 
of the nucleated bubbles will be analyzed in Section 3. In Section 4 we shall discuss 
the implications of the recycling universe model for the question of the beginning of the 
universe. The superlarge-scale structure of a recycling universe will be studied in Section 
5 using the methods of stochastic inflation. The issue of predictions will be discussed in 
Section 6, and our conclusions will be summarized in Section 7. 



II. MODELS 

In all realistic inflationary models, the potential of the inflaton field <^ is required to 
have a sufficiently flat slow-roll region in which 

\V"{(P)\<^H^ (3) 

Here, H is the expansion rate and is the horizon size corresponding to the vacuum 
energy V{(f)), 

H'^ = 87r1/(0)/3. (4) 

On the other hand, Coleman-de Luccia-type solutions for vacuum bubbles exist only when 
the potential is sufficiently curved near the barrier separating true and false vacua ||29|| , 

\V"{<P)\>H^ (5) 

The meaning of this condition is easy to understand. The bubble wall thickness is 5 ~ 
I and if (§) is not satisfied, then the wall is much thicker than the horizon. Such 

walls cannot exist as coherent structures and are spread by the expansion of the universe. 
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False vacuum bubbles of Coleman-de Luccia type are, therefore, impossible if the slow 
roll condition is valid everywhere in the inflationary range of 0. 

A similar problem arises in the "open" inflation scenario, where false vacuum decay 
through bubble nucleation is followed by a period of slow roll in bubble interiors @J^. 
One way to deal with this problem is to consider a two-field model, with one field doing 
the tunneling and the other doing the slow roll 0. The potential can be chosen as 

U{xA)=V,ix) + x'V2i<P). (6) 

Here, x is the tunneling field and the potential Vi{x) has the form as in Fig. 2 with a 
metastable minimum at x = and a true minimum at x = Vx- '^^^ ^^^^ potential U{x, <f)) 
is independent of (p (has a flat direction) at x = 0, and as a result the expansion rate 
in the false vacuum is also independent of 0. The potential V2(0) is assumed to have 
a slow- roll range and a minimum at (p = t]^ with V2(?7(/,) = 0, at which thermalization 
eventually occurs. In a variant of this model P], the two fields can be taken to represent 
the radial and angular parts of a single complex field, $ = xe*"^. 

The stage for open inflation is set by inflating false vacuum with energy density 
Pf ~ ^i(O)- Nucleating bubbles expand into this background, but because of the high 
expansion rate of the intervening false vacuum regions, bubble collisions are rare. The 
interior geometry of each bubble is that of an open Robertson- Walker universe. The 
bubbles have different initial values of (p, and if this value falls in the slow-roll range of 
the potential, then there is a period of inflation inside the corresponding bubble. Inflation 
is followed by thermalization and standard cosmological evolution, but since we assumed 
a non-zero cosmological constant, the bubble interiors are eventually dominated by the 
true vacuum energy, pt = Vi{ri^). False vacuum bubbles will now be formed in the true 
vacuum background, resulting in the endless succession of stages of the recycling universe. 

It should be noted that tunneling back to the false vacuum can occur not only from 
the true vacuum {x,4>) = {flxi'l<t>)i but also from the slow-roll, as well as radiation and 
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matter-dominated periods. In fact, the rate of false vacuum bubble nucleation is expected 
to be the highest during the slow-roll inflation [due to the higher energy density at that 
time, see Eq.(|^)]. However, since all these periods last only for a finite time, and the rate 
of false vacuum bubble formation is extremely low, only a tiny fraction of the co-moving 
volume will be affected by such processes. On the other hand, the true-vacuum-dominated 
stage persists indefinitely, and practically all the co-moving volume will be recycled by 
bubbles nucleating in the true vacuum. 

An alternative to the two-field model is a model of a single scalar field (f) with a 
potential of the form shown in Fig.3 ||^. The false vacuum at = is separated from the 
slow-roll region by a sharp barrier. (The coexistence of flat and highly curved regions in 
the same potential is a somewhat unnatural feature of this model). The field (j) tunnels 
through the barrier and after a period of slow roll, ends up in the true vacuum, which we 
assume to have a non-zero energy density, pt > 0. An important difference of this model 
from that of Eq. (P) is that now, in order to tunnel from true to false vacuum, the field (j) 
has to go across the whole slow-roll region. False vacuum bubbles will therefore consist of 
a false vacuum core surrounded by a domain wall, which is in turn surrounded by layers 
of slow roll, radiation, and matter-dominated regions (see Fig. 4). 

Intuitively, we would expect that the nucleation of such a complicated structure 
should be extremely unlikely, and thus the nucleation rate of false vacuum bubbles in 
this one-field model should be strongly suppressed compared to the two-field model (H). 
However, according to Eq. (|^) this would not be so. The bounce solution depends only on 
the form of the potential V{(f)) in the vicinity of the barrier. Thus the tunneling rate given 
by the expression (|^) is essentially independent of what happens in the slow-roll region 
and would not be much different if that region was absent. It should be remembered 
though that the use of Euclidean methods in de Sitter space has never been justified from 
first principles, and therefore the results obtained using these methods should be taken 
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with caution |8HlO|. This issue needs further investigation, but we will not attempt to 



address it in the present paper. The specific value of false vacuum bubble nucleation rate 
will not be important for our conclusions, as long as this rate is non-zero. 

If the barrier in the inflaton potential is too wide to satisfy the condition (j^) , then 
the Coleman-de Luccia instanton does not exist. However, there is always a homogeneous 
Hawking-Moss instanton [0] in which takes the value (pb corresponding to the top of the 
barrier. This instanton is usually interpreted as describing quantum tunneling from false 
vacuum to the top of the barrier in a horizon-size region. (Coleman-de Luccia instanton 
reduces to that of Hawking and Moss as the barrier width is increased). In the spirit of 
Lee and Weinberg |TT|, we can interpret the same instanton as describing tunneling from 



true vacuum to the top of the barrier. The corresponding nucleation rate is 

r (xexp[-SHM + St], (7) 

where Shm = —3/8pb and pb = V{(f)b)- 

In models of "new" inflation, the generic potential is illustrated in Fig.5. There is 
no barrier in this case, but still there is a Hawking-Moss instanton with (pb corresponding 
to the maximum of the potential. There are also approximate homogeneous instanton 
solutions with sufficiently close to the maximum. Such approximate instantons also exist 
in models of "chaotic" inflation where the potential may have no maxima. A constant 
field is a good approximation as long as the evolution of is slow on the Hubble scale 
H-^, that is, in the slow roll range. The Euclidean expression for the tunneling rate from 
true vacuum to in this range is given by Eq.(|^ with Shm = — 3/8V(0). 

Once again, we find the formula somewhat suspicious, and emphasize the need 
for a derivation of the nucleation rate without relying on Euclidean methods. Linde |]10 



has given an estimate of the probability of tunneling to false vacuum on rather general 
grounds. In de Sitter space, a field fiuctuates on scales bigger or comparable to the 
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Hubble radius around a local minimum (pt of the potential with amplitude given by 



CT^ =< id) - (hf)'^ >^ 



Strictly speaking, the result is only true for a free field and for small fluctuations, so 
that the mass m and the expansion rate H are well defined. Extrapolating to the case 
of an interacting field and ignoring the gravitational backreaction of the fluctuations on 
the expansion rate, the probability that a region of size ~ H^^ will tunnel from 0t to a 
different value sufficiently close to the maximum of V{(j)) can be estimated as 

Pocexp[-(0-0i)V2a2]. 

For the case of a quadratic potential with m? « and — V{(j)t) ^ the 
exponent in the previous expression reproduces the exponent [—Shm + St\ which appears 
in Eqn. ©, with Shm = -3/8^(0) |T0|. 

In any case, we expect the rate to be non-zero in the general case. The reason is 
simply that the nucleation is not forbidden by any conservation laws, and thus should 
have a non-zero probability. 



III. FALSE VACUUM BUBBLES 

To study the geometry of false vacuum bubbles, we shall first assume that the bubble 
wall thickness is much smaller than all other relevant dimensions of the problem. The 
wall can then be treated as infinitely thin, and the spacetime regions on the two sides 
of the wall are de Sitter spaces of different vacuum energy. We now briefly review some 
properties of de Sitter space. 

It is well known that de Sitter space can be pictured as a hyperboloid embedded in 
a flat 5-dimensional spacetime, 
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where ( = {(^, C^, is a 3- vector. A section of the hyperboloid by the wr-plane is shown 
in Fig. 6. The Euchdeanized de Sitter space, which is used for constructing instantons, is 
obtained by analytic continuation r = iTE, 

+ + r| = H-\ (9) 

Geometrically, this is a 4-sphere of radius H~^. 

It will be convenient to use the Robertson- Walker flat coordinates in which the de 
Sitter metric takes its most familiar form, 

ds^ = dt^-e^'''dx.\ (10) 

These coordinates are related to the hyperboloid coordinates by 

r = H-^ sinh(m) + ^Hx^e^\ (11) 

w = H-^ cosh(m) - ]-H^^e^\ (12) 

C = xe^*, (13) 

which can be inverted to give 

t = E-^\^\E{w^r% x=— p^. (14) 

Constant-t surfaces are obtained as intersections of the hyperboloid with null hyperplanes 
w + r = const (see Fig. 6). The surface t = — oo corresponds to w + t = 0, and thus the 
coordinate system (|TD|) covers only half of de Sitter space. 

Let Hj^ and Hf^ be de Sitter horizons corresponding respectively to the false and 
true vacuum energy densities, pj and pt- Clearly, pt < pf and Hj^ < H^^. The thin wall 
approximation requires that the wall thickness be much smaller than Hj^. In this case, 
the Coleman-de Luccia instanton for bubble nucleation can be obtained by matching two 
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4-spheres of radii Hj ^ and ^ (see Fig. 7). The two spheres are joint at a 3-sphere which 
represents the Euchdeanized worldsheet of the domain wall. Its radius Rq is determined 



by p/, pf, and the wall tension o |2T|. The 5-dimensional coordinates can always be chosen 



so that this worldsheet lies in a plane of constant w. In the figure it is = — -Dt, where 

D, = {H^' - Rlfl\ (15) 

The Lorentzian evolution of the bubble is given by two hyperboloids similarly matched 
along a constant-w plane (Fig. 8). Descriptions of both true and false vacuum bubbles can 
be obtained with an appropriate slicing of this spacetime by equal-time surfaces. 

In the case of false vacuum bubbles, equal-time surfaces can be chosen to be the 
surfaces w + t = const (Fig. 8). Then, each constant-t slice consists of a spherical region 
of false vacuum embedded in an infinite, spatially-flat region of true vacuum. Since the 
spatial geometry of these slices is flat both inside and outside the bubble, the volume that 
is removed from true vacuum by the appearance of the bubble is equal to the volume 
of false vacuum which replaces it (this would not be the case if we used closed spatial 
sections, for instance). We can use the coordinates ( p!0D with H = Ht to describe the 
exterior true- vacuum region. The wall worldsheet is at w = — -Dj, and from Eq.(|12]) the 
radius of the bubble at time t is 

R}{t) = H-^e^""'' + 2Ai^te^** + 1), (16) 

where R(t) = |x(t)|e'^**. We see that the radius approaches the horizon size Hj~^ at 
t — s> — oo. The bubble wall accelerates in the direction of the false vacuum, so that its 
co-moving radius \^(t) \ is contracting, but the physical radius grows exponentially due to 
the expansion of the universe. A conformal diagram for the bubble spacetime is shown in 
Fig.9. 

Since R{t) is a monotonically growing function of time, there is no "bounce" moment 
at which one can say that the nucleation occurs. The situation here is similar to that for 
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nucleation of topological defects in de Sitter space and as in the latter case, we shall 



regard the bubble (|T6|) as "formed" at t ~ 0, when its radius begins to grow exponentially. 



For the trajectory given in Eqn. (16), the region of true vacuum which has been 



removed and replaced by true vacuum was centered at the point (t = 0, x = 0) (see Figs. 
6 and 8). We shall refer to this point as the center of symmetry of the bubble trajectory. 
(Notice that this point may not belong to the actual classical spacetime, because it is 
precisely in the region where bubble nucleation takes place, but it does belong to the five- 
dimensional embedding space.) Performing Lorentz transformations in the embedding 
space, we can obtain bubbles whose center of symmetry is at any time t = to and at any 
location. This gives 

Rj(t) = i7-2[e2J^t(i-io) + 2DtHte^'^'-'''^ + 1]. (17) 

We can think of these as bubbles "formed" at t ~ to- The asymptotic behavior of the 
bubble radius at large times is 

Rfit) ^ H~^e"'^'-'°\ t-to^H-\ (18) 

To describe a true vacuum bubble in a false vacuum background, we choose equal- 
time surfaces to be w — r = const. Using the coordinates (|1^) with H = Hj to describe 
the exterior of the bubble, we find that the bubble radius at time t is given by 

Rf{t) = Hf{e^"f' - 2DfHfe"f' + 1), (19) 

where 

Dj = iHj'-RlY/'. (20) 

In contrast to the false- vacuum bubble case, the radius ( |T9D has a minimum, Rmin = Ro, 
at tn = Hj^ In(DfHf), and we can regard this time as the moment of bubble nucleation. 
Eq. (|T9|) can then be rewritten as 
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R^tit) = /)J[e2^^/(*-*») _ 2e^/(*-*")] + Hj^, (21) 
with the late-time behavior 

Rt{t) ^ Dfe"f'-'-'"\ t-tn^ Hj\ (22) 

Note, however that 

Rt{t) ^ Hj^e^f'^'-'°\ t - to > Hj\ (23) 

where we denote by to the location of the center of symmetry of the bubble wall as seen 
from the outside [see the discussion around equation (0)]. 

Our main interest in this paper is in models with pt ^ P/ • In such models, the radius 
of false vacuum bubbles is R{t) > H^^ ^ -^7^' thus the thin wall approximation can 
be used to describe the bubble evolution even when the wall thickness is 5 ~ ^J^- (Note 
however that in this case the thin wall approximation breaks down for the instanton itself 
and for the early evolution of true vacuum bubbles). 

In a single-field model of open inflation, with a potential as in Fig. 3, the conformal 
diagram for a false vacuum bubble is shown in Fig. 10. 

So far in this Section we assumed that the true vacuum has a positive energy density. 
If the vacuum energy is in fact zero, then the horizon radius in thermalized regions keeps 
growing with time, and false vacuum bubbles eventually come within the horizon. The 
bubbles are then seen as black holes from the outside. The bubble nucleation in this 
case is similar to the quantum creation of baby universes, as discussed in Refs. [Q,^,|T^ . 
Black holes eventually evaporate and baby universes pinch off. 

IV. DID THE UNIVERSE HAVE A BEGINNING? 



Assuming that some rather general conditions are met, it was shown in Ref. ||T5 
that inflationary models cannot be geodesically complete to the past, that is, they require 
some sort of a beginning. The assumptions that lead to this result are the following. 
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A. The universe is causally simple ||36|| . 

B. The universe is open. 

C. The null convergence condition [^]. 

D. The finite past-volume difference condition. 

The first two of these assumptions do not appear to be crucial for the proof, and 



extensions of the theorem have been obtained to some closed universes [161 and to some 



universes with a more complicated causal structure [17]. 

The null convergence condition is closely related to the weak energy condition, which 
requires that the energy density is non-negative when measured by any observer. Clas- 
sically, this is satisfied by all known forms of matter, including a relativistic scalar field, 
but violations of the null convergence condition are possible as a result of quantum fiuc- 
tuations. Such violations tend to occur in the infiating regions of spacetime whenever 



quantum fiuctuations result in a local increase of the expansion rate, dH/dt > [|4 



They are sufficient to invalidate the theorem in models characterized by a substantial 
variation of in the range of where quantum fiuctuations are non-negligible. This 
includes all models of "chaotic" infiation, but not some open and "new" infiationary 
models. 

Turning now to the effects of recycling, we shall argue that they invalidate both 
assumptions C and D, so that the theorem as it stands does not apply to any inflationary 
model. In the case of assumption C, the reason is the same as before: quantum fluctuations 
from thermalized regions back to the inflating phase increase the expansion rate in the 
affected regions of space, and the null convergence condition is violated. 

The flnite past-volume difference condition {D) can be formulated as follows. Given 
a point P in the inflating region and a point Q to the past of P, consider the difference 
of their pasts. This is a spacetime region including all points to the past of P, but not of 
Q. The condition D requires that the spacetime volume of this region should be flnite. 
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The original motivation for this condition was based on the picture of eternally inflating 
universe without recycling, as illustrated in Fig.l. It can be shown that thermalization 
surfaces, which separate inflating and thermalized regions of spacetime, are spacelike 



surfaces [15]. Therefore, if Q is a point in an inflating region, then, disregarding recycling, 
all points in its past must also lie in the inflating region. For inflation to persist from Q 
to P, no thermalized regions should be formed in the difference of the pasts of the two 
points. Now, it seems plausible that there is a zero probability for no thermalized regions 
to form in an infinite spacetime volume. Then it follows that condition D is necessary 
for inflation to persist in the future time direction. This condition is difficult to satisfy, 
since the spacetime region representing the difference of the pasts of the two points has 
an infinite extent along the past-directed null geodesies. 

In a recycling universe, this logic does not apply, since points in inflating regions 
can have thermalized regions in their past. In this case, the spacetime region relevant 
for persistence of inflation between a pair of points is not the entire difference of their 
pasts, but only the part of this difference which lies in the same inflating region as the two 
points. In other words, it is the part of the difference of the pasts which is to the future of 
the "nucleation surface" (see Fig.9). The volume of this region is obviously finite. Hence, 
Assumption D is not suitable for a recycling universe. 

In models with a vanishing true vacuum energy, only a small fraction of thermalized 
volume gets recycled. However, there still appears to be a possibility that the universe has 
a nested structure, with all inflating regions originating as quantum fluctuations inside 
thermalized regions. 

We thus conclude that the theorems of Refs. [|15H T7| no longer apply when the recy- 
cling nature of the universe is taken into account. This may open the door for constructing 
non-singular, steady-state inflationary models. We emphasize, however, that our analysis 
does not imply that such models do in fact exist. It has been argued in Refs. [^,0] that 
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inflation, when continued to the past, is necessarily preceeded by a period of contraction, 
as in the exact de Sitter space. During this period, the thermahzed regions would merge, 
the density perturbations would grow very fast, and the universe would rapidly reach a 
grossly inhomogeneous state from which it is not likely to recover. The arguments in 
3^,34] do not rely on weak energy or finite past-volume difference conditions, and may 



possibly be extended to the case of a recycling universe. These arguments, although 
suggestive, fall short of a proof, and the problem requires further investigation. 



V. STOCHASTIC FORMALISM 

A quantitative description of the recycling universe can be given using the formalism 
of stochastic inflation developed in Refs. [pi|8|,pj|. A straightforward extension of this 
formalism will be required, and to introduce the necessary modiflcations, we shall flrst 
consider the Lee- Weinberg model [|ll[] with a potential as in Fig.2. 



A. Lee- Weinberg model 

Consider an ensemble of co-moving observers whose world lines are orthogonal to 
some spacelike hypersurface S. Let P/(t) and Pt{T) be the fractions of observers in false 
and true vacuum, respectively, 

Pf{T)+Pt{T) = l, (24) 

where r is the proper time along the observer's world lines measured from their intersection 
with S. The time evolution oi pf and pt is described by the system of equations 

dpf/dr = -KfPf + Ktpt, (25) 



dpt/dr = -Ktpt + KfPf. 
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(26) 



Here, Kf (k^) is the probability, per unit time, for an observer who is presently in the false 
(true) vacuum to find himself within a true (false) vacuum bubble. 

^From Eq. (pl|), we see that a false- vacuum observer will be affected only by bubbles 
nucleating within a sphere of radius Df centered on that observer. The bubbles take time 
~ Hj^ to traverse this distance, but in the stochastic inflation formalism we shall be 
interested in quantities smeared over a spacetime scale ~ H^^, so we shall disregard this 
time delay and write 



nf^Tf^fD% (27) 



where Fj"^ is the rate of bubble nucleation per unit spacetime volume. 

In an expanding universe, however, the rate of nucleation per unit spacetime volume 
has to be defined with some care. This is because this rate depends on what we choose 
as the nucleation time, which is not always sharply defined (especially when the size of 
bubbles becomes comparable to the Hubble radius.) For instance, we can change our 
definition of nucleation time from t„ to to, where t„ is defined by Eqn (^) as the time at 
which the physical radius of the bubble reaches its minimum value i?o, and to is defined as 
the center of symmetry of the bubble wall trajectory as seen from the outside (as explained 
in Sect. III). In changing the definition of nucleation time, we must simultaneously 
change the definition of nucleation rate per unit voulme, because the physical volume has 
increased by the amount exp[3 Hf {to — tn)] in the intervening time. Distinguishing by their 
superindex the rates associated with both choices of nucleation time, we have 

rf^ =Tfexp[3Hf{to-t,,)]. 

Hence, Eqn. ( pTj) can be rewritten as 

K,«rf|ff-. (28) 

In the case of false vacuum bubbles the physical radius grows monotonically with 
time, so the analogue of t„ does not exist. We shall adopt the convention identifying the 
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nucleation time with to in Eq. (|T3). Then the co- moving region affected by the bubble is 
a sphere of radius H^^, and we can write 

K,«rf'|ff-, (29) 

where vf'^ is the corresponding nucleation rate. We note that the radius of the affected 
region and the rate would both be modified with a different choice of nucleation time, 
while Kt would remain unchanged. 
The solution of Eqs.(p5|),(p^) is 

Pf{T)=pf + Ae-'^, (30) 

Mr)=pf-Ae-'\ (31) 

where the constant A is determined by the initial conditions, b = kj + Kt, and p^p , pf* is 
a stationary distribution defined by 

pT+pf' = l- (33) 
We see that, regardless of the initial conditions, the probability distribution rapidly ap- 



proaches the stationary distribution (|3^). 



Let us compare the distribution (|32D with that discussed by Lee and Weinberg 



TT|. The distributions agree provided that the nucleation rates which appear in their 

(0) 

t/r 



expressions are taken as Fwl. Lee and Weiberg also argue that the ratio (R^) can be given 



to one loop order as 



- = S^^P[-(3/8)(pr^-p7^)], (34) 

Kf Ut 

The absence of determinantal prefactors in this expression is justified by the fact that 
the bounce solution for true and false vacuum decay are the same, hence the primed 
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determinants corresponding to fluctuations around the bounce cancel out in the ratio 
of rates. The effect of determinants corresponding to fluctuations around the true and 
false vacuum background instantons is just to renormalize the values of the corresponding 
effective cosmological constants. Hence, in (|3^) , the vacuum densities pt/f should be 
understood as their 'one loop corrected' values. 

In Ref. the nucleation rates for true and false vacuum bubbles were studied 
in the case when the gravitational backreaction of the bubble is ignored, so that the 
background geometry was taken to be an exact de Sitter space with the same Hubble 
constant H inside and outside te bubbles. It was found that the number distribution 
of true or false (t//) vacuum bubbles centered around the point (to,x) could be written 
as dMt/f = \Xt/f\exp (SHto) rf^x dtQ. Here Xt/f = Aexp {~Bt/f), where A is a primed 
determinant which is the same for true and false vacuum bubbles, and Bt/f = Sb — St/f 
is the difference between the bounce action and the background Euclidean action. The 
expression for dM is proportional to the physical volume element at time to? given by 

t/f \ with the rates vfj^^ 



exp (SHto) (i^x, so we can identify |At//| with the rates F^l defined above. Therefore, we 



have {t[^^ /vf^) = exp [St — Sf], in agreement with Eqs. (|32| , |34D . 



We introduced j9/ and pt as fractions of co- moving observers in false and true vacuum. 



respectively. When the stationary distribution ([3^ ) is reached, an alternative interpreta- 
tion will also be useful. The world line of each observer will repeatedly cross between true 
and false vacuum regions, and the quantity pj*"* {pf'^) gives the fraction of proper time 
the observer spends in false (true) vacuum. 

Instead of the proper time r, one can use some other time variable, t, along the 
observer's worldlines. A possible choice is 

dt = H'^{T)dT, (35) 

with a = corresponding to the proper time and a = 1 to the "scale factor time". For 
an arbitrary a, the evolution equations still have the form (125|) , (^) , with r replaced by 
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t and 




-a-3 



(36) 




-a-3 



(37) 



The stationary solution now is 



Pf 



Pt 



( 




(38) 



The a-dependence of (^) can be easily understood: pf {pt) is proportional to the amount 
of time spent by a co-moving observer in false (true) vacuum, and if the time variable is 
changed as in (|35D, the ratio pj/pt is modified by a factor (Hf/Ht)'^. 



Let us now consider a two-field model of the type (P), except we shall assume that 
the false vacuum at x = corresponds to a single point, rather than a fiat direction, in 
the field space. This is the case, for example, in models where x ^^id represent the 
radial and angular parts of a single complex field, $ = xe*"^. We shall assume further that 
the effective potential for the field is of the "new" inflation type and has a slow roll 
region ^l^"* < < ^l^"* . Finally, to simplify the discussion, we shall disregard the evolution 
between the end of slow roll and true vacuum domination. That is, we shall assume that 
when the field rolls down to (f)i^\ it gets directly to the true minimum of the potential 
with energy density where j = 1, 2. We shall refer to 01"' ■* as "thermalization points" 
and to the corresponding minima of the potential as the first and the second true vacua, 
TV{1) and TV{2). (If p[^^ ^ pf\ then of course only one of these vacua is a truly true 
vacuum). The false vacuum will be abbreviated as FV . 

Once again we introduce an ensemble of co-moving observers and define Pti{t), Pt2{t) 
and Pf{t) to be the fractions of the observers in TV{1), TV{2) and FV, respectively. We 



B. A more realistic model 
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also define V{(f),t)d(j) as the fraction of observers wlio are, at time t, located in slow roll 
regions with the infiaton field between (j) and (f) + d(f). We can now combine the analysis 
in the preceeding subsection with the standard formalism of stochastic inflation to obtain 



the system of equations describing the evolution of our model: 

dV dJ 

^^ = "90"'^^^'^^^ + ''^^'^^^^' ^^^^ 

= -i^nPti + t/iP/ - A, (40) 

= -fit2Pt2 + fif2Pf + J2, (41) 

^ = -(/t/1 + K,f2)Pf + K,tlPtl + K,t2Pt2 -Pf I K,fi(l))d(f) + [ K f {(t))V d(f) . (42) 



Here, Kfj corresponds to tunneling from FV to TV{j), Htj to tunneling from TV{j) to 
FV; they are given by Eqs.(PB|), ( P7[ ) with an extra index j added to the appropriate 
quantities. K,f{(f))d(j) corresponds to tunneling from FV to a value (p in the interval d(j) in 
the slow roll region, and iif{4>) to tunneling from a slow roll region with a given value of 
to FV. By analogy with (|36D, (|37|) we can write 

'^/(0) = ^f(0)fi/7°-^ (43) 

^/(0) = ff(0)yi/-"-^(0), (44) 

To simplify the equations, we have disregarded tunneling between TV{1) and TV{2) and 
between TV{j) and the slow roll region. These effects can be easily included if necessary. 
The current J{4>,t) in Eq. ( pQ] ) is given by 

J(0,t) = -Di-^(0)|-[D^(0)P(0,t)] + i;(0)W,t), (45) 
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where the first term on the right-hand side describes the "diffusion" of the field due to 
quantum fluctuations, with a diffusion coefficient 

D{<j>) = H'-^{<j>)/87r', (46) 

the second term describes the classical "drift" with velocity 

v{(P) = -H-^{4>)H'{^)/4n, (47) 

and 

H\(f)) = 87r1/(0)/3. (48) 

The integration in Eq.(|4^) is from ^1^^ to (f)^*\ and the quantities Jj in Eqs.(^),(Pl) are 
defined as Jj{t) = J{(f)i^\t). The normalization condition 

j Vd(t) + Pf+pti+Pt2 = l (49) 
is preserved by the evolution equations (|39D-(^2D. 



The parameter (3 in Eq. ( ^5]) for the current represents the factor-ordering ambiguity 
in the diffusion equation with a position-dependent diffusion coefficient. The choices 
(3 = 1/2 and (3=1 are usually referred to as Stratonovich and Ito factor ordering, 
respectively. 

The boundary conditions for Eq.(P^ are 

^[D^(0)P(0,t)]^^^o)=O. (50) 

They ensure that, once the field rolls down to 4)"J\ it does not diffuse back to the 
slow-roll regime, but rather stays in the true vacuum TV{j) until it tunnels to FV . 

The system of equations (^)-(|4^) can be written symbolically in the operator form, 

^=MP. (51) 
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where the "vector" P(t) is P = {pf,Pti,Pt2,'P{4>)}- With an appropriate discretization of 
(j), this can be rewritten in the form of a "master equation" , 

r/P 

— ^ = EK-^. - ^..^.) = E M^,Py (52) 
i j 

Each quantity is positive and has the meaning of the transition rate from state j to 
state i. The matrix Mij can be represented as 



Mij = Wij - 6ij Wki (53) 

k 



and has the properties 



M,, >0 {i^j), (54) 
^M,, = 0. (55) 

i 

The latter property ensures the conservation of probabihty, X^i-Pj = 0. It also indicates 
that the matrix M has a left eigenvector Q = {1, 1, 1, ...} with zero eigenvalue, QM = 0. 
Since M and its transpose have the same eigenvalues, it follows that M should also have 
a right zero eigenvector, 

MPo = 0, (56) 

indicating that our system of equations has a stationary solution. 

The familiar method of solving Eq.(|5^ using a decomposition in eigenvectors can- 
not, in general, be applied because the matrix M is not generally symmetric. However, 
some properties of its solutions can be derived using only Eqs. ([5^) , (^) . We shall assume 
that the matrix M is irreducible (otherwise, the master equation ( ^I]) describes several 
independent processes which can each be described by a master equation with an irre- 
ducible matrix M). Then it can be shown that (i) the zero eigenvalue, 70 = 0, is 
non-degenerate, (ii) all components of the corresponding eigenvector Pq are non-negative, 
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(iii) all other eigenvalues of M satisfy Re7„ < 0, and (iv) the asymptotic behavior of P(t) 
at late times is 

P(t ^ oo) = Po. (57) 

In other words, the stationary solution is unique, and all solutions approach this stationary 
solution at t ^ oo. Although these results have been rigorously derived only for a finite 
set of Pn, we shall assume that they are still valid in the continuum limit. 

If all eigenvalues of M are non-degenerate, then M can be diagonalized, and the 
general solution of (0) can be written as 

oo 

P(t) = Y: Pne^-\ (58) 

n=0 

where P„ are eigenvectors of the operator M, and 7^ are the corresponding eigenvalues, 

MP„ = 7„P„. (59) 

Since M is real, its eigenvalues and eigenvectors come in complex conjugate pairs. In the 
case of degenerate eigenvalues, the solution is more complicated ^ . 



In the absence of recycling, Ktj = i^f{<P) = 0, and the stationary solution of the sys- 
tem (PP|)-(P2|) is trivial: ptj = const, pf = V{4>) = 0. The standard analysis of stochastic 



inflation P,P|,P^ has been done for a slow-roll inflation without a metastable false vacuum. 
Then Eq.(^) reduces to a Fokker-Planck equation for V{(f),t), 



It can be shown ||41|] that, with an appropriate choice of a scalar product, the differential 
operator Ai is Hermitian. Hence, all its eigenvalues are real and the eigenvectors form 
a complete orthonormal set. An eigenvector expansion of the form (|58D is then always 



possible, and the asymptotic behavior of V{t) is 

P(0,t-^oo) = /(0)exp(7„t). (61) 
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Here, 7^ < is the largest eigenvalue of A^. In this standard approach the distribution 
V{(f), t) is not normalized: the observers who left the slow roll range through the boundaries 
at ^l"'^ never return, and / Vd(f) decreases exponentially with time. 



C. The fractal dimension 



In the case of "new" inflation without recycling, the inflating part of the universe 
represents a self-similar fractal of dimension [| 



d = 3-^^H^-\ (62) 

Here, as in Eq.(^), 7^ is the largest eigenvalue of the Fokker-Planck operator JH, and 
Hm is the expansion rate at the maximum of V{(f)) |^2|. For a co- moving sphere of radius 
R centered on a point in the inflating region, the inflating volume within the sphere is 
V(-R) oc R^, which is a fraction 

f{R) oc i?'^-^ (63) 

of the total volume of the sphere. As the sphere expands, R oc exp{H^"'t) , this fraction 
decreases as / oc exp(7mt), and vanishes at t ^ 00. Hence, the inflating region represents 
a set of measure zero in the limit t ^ 00. 

On the other hand, the inflating part of the volume in a recycling universe is con- 
stantly replenished by tunneling from the true vacuum. As a result, the inflating region 
occupies a non-vanishing fraction of the total volume, so that V{R) oc i?^ and d = 3. 

We note, however, that a recycling universe does contain fractal regions of dimension 
d < 3. Take for example the Lee- Weinberg model of Section A and consider a co-moving 
volume which is initially filled with FV. What remains of this FV in the limit t 00 is 
a fractal of dimension 

d, = 3-TffHj\ (64) 
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All the remaining part of the volume is occupied by true vacuum bubbles, but what 
remains of the TV inside the bubbles at t — > cxd is also a fractal of dimension 

d, = 3-Tr^fH-\ (65) 

The FV bubbles inside each TV bubble have dimension df, and they are in turn filled by 
TV bubbles of dimension df. The fractal structure of realistic models is of course more 
complicated. 



D. Choosing the factor ordering 

One of the problems with interpreting the results of the stochasic inflation formalism 
is the dependence of these results on the choice of the time variable t and on the factor 
ordering in the Fokker-Planck equation (|39|) , (45) |19;27]. We have parametrized these 



choices by the parameters a and (3. Now we are going to argue that there is a preferred 
choice of f3 which allows at least a partial resolution of these problems. 

As we discussed in Section IV. A, the stationary distribution Pq gives the fraction 
of time spent by a co-moving observer in false and true vacua and in different parts of 
the slow-roll range. This distribution should of course depend on how we define the time 
variable, but the dependence should be trivial, 

P(0)(xi/"(0), Pf^Hf, Ptj^H^y (66) 

To ensure that this is indeed the case, we can require that in the stationary version of 
Eq.(§T]), 

MP = 0, (67) 

all components of P appear only in combinations iJ~"(0)7^(0), Hj"pf, H^°'ptj, and that 
there is no other dependence on a. This fixes 
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13 = 1, (68) 

which corresponds to Ito factor ordering. 

It should be noted that the family of factor ordering choices parametrized by j3 in 
Eg. (^51) does not include all possibilities. Although our requirement (|66D is sufficient to 



fix j3 uniquely, it does not determine a unique factor ordering in the general case. For 
example, we could replace the diffusion term in (HSl) by 



The condition ( |BBD is satisfied for an arbitrary, a-independent function /i(0). 



E. Some solutions 

The system of equations (^)- (|4^) looks rather intimidating, but stationary solutions 
of this system can actually be found in some special cases. 

Let us first assume that tunneling from slow roll to false vacuum can be neglected, 
that is, K/(0) = 0. This means that the evolution proceeds along the path FV slow 
roll TV —>■ FV — > possibly with occasional tunneling directly from FV to TV. 
Then, with /? = 1, the stationary version of ( PUj) can be written as 

9^J(0) = K{^)pf, (70) 

where 

J{<f>) = -d4D{cl>)Vm + v{<jy)V{<jy). (71) 
This is easily integrated to give 

(72) 



P(0) 
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where we have used the expressions (^),(^7|) for D{(f)) and f (0). 

The integration constants Ci, C2 can now be found from the boundary conditions 
(0) (with f] = 1). Thus we obtain the distribution V{(f)) in terms of pf. Evaluating 
Jj = J{(t)^P) and substituting in Eqs. (^0D , (|4l|) , we find ptj in terms of p/. Finally, p/ is 



found from the normalization condition (49). [Note that pj cannot be found from Eq.( 
which is a linear combination of the preceeding three equations (^)-(0)]. The resulting 
expressions are rather cumbersome and we shall not reproduce them here. 



As another example, we take a potential V{(j)) of the form considered in Ref. [21 



which consists of a fiat portion where -ff(0) = const, Hf{4>) = const, Kf{4>) = const, sur- 
rounded by two regions with a relatively large slope where the diffusion term is negligible. 
In the flat range of (p, the Fokker-Planck equation is trivially solved. In the difFusionless 
regions. 



(73) 



and a straightforward integration gives 



P(0) 



(1) 



where 



(1) V{(f)') 



(74) 



(75) 



for the region bounded by ^l^"*, and similarly for the second region bounded by 
The integration constants and the values of p/ and ptj can be determined by matching 
the solutions at the boundaries between the fiat and diffusionless regions and by using 



Eqs.(^),(0) and the normalization condition (^ 



In the general case, the solution (|7^ should still apply in the range of (f) sufficiently 
close to the thermalization points, where diffusion is negligible. If tunneling between this 
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range and the false vacuum is unimportant, then the solution takes a particularly simple 
form. 



(76) 



(2) 

and similarly for the range of near 0* . The constants Ci and C2 can be determined 
only after solving the equation in the entire range of 0. 

The distribution (|76D in the diffusionless range of should be compared with the 



corresponding distribution in the absence of recyclings ||2T 



W) = C^^exp 



H'{ct>) 



H' 



(77) 



This can be drastically different from (|76D even if the tunneling probabilities are very small. 
Hence, there is no continuous transition between recycling and no-recycling regimes in the 
limit of vanishing tunneling probabilities. 



F. Physical volume distribution 

The function P(t) characterizes the distribution of co-moving volume between false 
and true vacua and different values of in the slow roll regime. One can introduce a 



similar function for the physical volume distribution [^3H45[|, 



P(t) = {V{M,pf{t),pn{t),p,2m- (78) 

It is defined so that P(0, t)d(j) is the physical volume occupied by slow roll regions with 
in the range d(f) at time t, etc. The distribution P{t) satisfies a modified master equation, 

//p _ _ 

— = MP + 3H^ "P = MP, (79) 
at 

where the operator H^^° is given by 

H^-" = dtag{H'-^{<p), H}-^, HI^, H],--}. (80) 
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In an infinite universe, the volume is of course infinite, but the distribution P(t) can be 
defined on a fixed co-moving part of the universe. The form of the distribution at large t 
is independent of the choice of the co-moving region. 

In the discretized version of Eg . (|79D , the matrix M does not have the property 
(^), and the standard theorems for the master equation do not apply. However, the 
following statements can still be proved using the Perron- Frobenius theorem about 



non-negative matrices: (i) M has a real eigenvalue 70 which is greater than the real parts 
of all other eigenvalues and which is bounded by 

3-ff^7n < 7o < S-ff^a", (81) 

where Hmax and Hmin are respectively the largest and the smallest values of if; (ii) the 
corresponding eigenvector Pq has non-negative components; (iii) 70 is non-degenerate if 
M is irreducible. The late-time asymptotic behavior of P(t) is 

P(t ^00)= Poe^"*, (82) 

In contrast to the co- moving distribution P(t), the physical volume distribution P(t) 
has a sensitive dependence on the choice of the time paramater a, which does not reduce 
to the trivial form (|6^) [1^,21,27]. The equation ( ffOD for P(t) is simplified if we choose 



the scale factor time, a = 1. In this case, " = 1, and the solutions of ([79|) and (|5l| ) 



are related by |19 



P„=i(t) =e^*P„=i(t). (83) 



G. Discussion 



The main conclusion of our analysis in this Section is that the distribution of co- 
moving observers in a recycling universe rapidly approaches a stationary form. This 
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asymptotic distribution can be obtained as the eigenvector of the "master" operator M 
with a zero eigenvalue, 70 = 0: 

MPo = 0. (84) 

The formahsm we developed here can be straightforwadly extended to include radia- 
tion and matter dominated periods between thermalization and true vacuum domination. 
One expects to find that the asymptotic distribution will still be stationary, with fixed 
fractions of co-moving volume occupied by radiation and matter-dominated regions. 

The picture in which co-moving "observers" move in an endless cycle between FV, 
slow roll, matter domination, and TV, may be oversimplified. It should be understood, of 
course, that no material observer is likely to survive the transition between TV and FV. 
Even if we think of an "observer" as an indestructible test particle, there seems to be no 
unique way to continue his world line into a nucleating bubble, since the surface at which 
we glue the bottom of the false vacuum bubble onto the true vacuum can be chosen in 
different ways. So we should probably think of our "observers" as being smeared over a 
horizon-size volume. 

Next, we note that density fluctuations produced during inflation (or generated by 
topological defects) result in the formation of bound objects during the matter-dominated 
era. Some of these objects collapse to black holes, and observers in matter-dominated 
regions have a flnite probability (per unit time) to end their world lines at black hole 
singularities. However, black holes eventually evaporate, giving back their volume to the 
true vacuum. Hence this effect would not alter our conclusions. The same happens with 
black holes that may spontaneously nucleate in false or true vacuum The rate of 



black hole pair production in true vacuum is proportional to exp[— l/8pt]. This rate is 
considerably larger than that for nucleation of a false vacuum bubble, which in the case 
Pt « Pf is proportional exp[— 3/8p(]. It is also possible that nucleated black holes may 



act as seeds for false vaccum bubble nucleation, as they do for true vacuum bubbles [|8 
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Finally, the inflaton potential V{(j)) can be of the "chaotic" inflation type, with the 
slow roll range of (p extending to Planckian energy densities. Then there is a finite prob- 
ability for an observer to get into this Planckian domain, where the classical description 
of spacetime breaks down. In the stochastic inflation formalism, this is accounted for by 
introducing a "Planck boundary" at some (j) = (j)p, such that V{(j)p) ~ 1. The loss of 
observers through the Planck boundary will generally result in 70 > and d < 3. 

The same phenomenon of loss of observers will also occur if some of the minima 
of the potential have vanishing or negative cosmo logical constant. Once some co-moving 
volume falls into one of these vacua it has no chance of being recycled. 



VI. PROBLEMS WITH PREDICTIONS 

Different thermalized regions of the universe are generally characterized by different 
values of the constants of Nature and of the cosmological parameters (such as the density 
parameter Q). In the model that we used as an example in Section V, the universe can 
thermalize into two types of vacua, TVl and TV2, and thus we have two possible sets of 
constants of Nature. The number of possibilities can, in principle, be much larger, and 
in some models the "constants" can even take values in a continuous range (examples 
are the effective gravitational constant in a Brans-Dicke-type theory and the density 
parameter in some models of "open" inflation |^) An intriguing question is whether or 
not we can "predict" which set of the constants we are most likely to observe. 

One can try to determine the probability distribution for the constants with the aid 
of the "principle of mediocrity" which asserts that we are "typical" among the civiliza- 



tions inhabiting the universe p0|j24| , |25||r8|| . Here, the universe is understood as the entire 
spacetime; our civilization is assumed to be typical among all civilizations, including those 
that no longer exist and those that will appear in the future. The probability for us to 
observe a given set of constants is then proportional to the total number of civilizations 
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in the corresponding type of thermalized regions. This number can be represented as the 
number of galaxies (which one can hope to estimate) times the number of civihzations 
per galaxy (which is left undetermined until the evolution of life and consciousness are 
better understood). Some of the constants, such as the cosmological constant or the den- 
sity parameter, are not expected to affect the chances for a civilization to develop in a 
given galaxy, so one can hope to determine the probability distribution for such constants 
without any biological input. 

In the case of a closed universe and finite (non-eternal) inflation, this prescription 
for calculating probabilities is unambiguous. If the universe is spatially infinite, one can 
simply use the prescription for a fixed (sufficiently large) co-moving volume. However, 
in an eternally inflating universe the spacetime volume and the number of civilizations 
are infinite, even for a region of a finite co-moving size. One can deal with this problem 
by simply introducing a time cutoff and counting only the number of civilizations Nj{tc) 
that appeared prior to some moment of time tc- Here, the index j refers to the type of 
thermalized region. The ratio of probabilities can then be defined as the limit 

Pi y 

— = hmt^^^——. 85 
P2 N2{Q 

One finds, however, that the resulting probability distribution is extremely sensitive to 
the choice of the time variable t |]49|JT9[| . This gauge-dependence casts doubt on any 
conclusions reached using this approach. 



An alternative procedure, suggested in |21|], is to introduce a cutoff at the time 
t^J\ when all but a small fraction e of the co-moving volume destined to thermalize into 
regions of type j has thermalized. The value of e is taken to be the same for all types of 
thermalized regions, but the corresponding cutoff times t[^^ are generally different. The 
limit e ^ is taken after calculating the probability distribution for the constants. It was 



shown in pT| , p7| that the resulting probabilities are essentially insensitive to the choice of 



time parametrization. However, the same problem appears in a different guise. Linde and 
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Mezhlumian have found a family of gauge-invariant cutoff procedures parametrized 
by a dimensionless parameter g, with q = corresponding to the e-procedure described 
above. This indicates that the invariance requirement alone is not sufficient to define the 
probabilities uniquely. Some additional requirements that can fix the parameter q have 
been discussed in . 

Now, recycling introduces one more difficulty. In the absence of recycling, co-moving 
regions could be uniquely characterized by the type of thermalized region they will evolve 
into. But in a recycling universe each co-moving region goes through an endless succession 
of different types of thermalization. Hence, the e-procedure cannot be implemented in its 
present form pU| . 

In the face of these difficulties, one could look for entirely different approaches 
to defining the probabilities. One possibility is to abandon the requirement of gauge- 
invariance and assert that there is a preferred choice of the time variable t. If this approach 
is taken, then there is, arguably, a good reason to take the scale-factor time, t = Ina, as 
the preferred choice [p|,|5l|. The only variables that can be used as clocks in an inflating 
universe are the infiaton field (j) and the scale factor a. The main requirement for a clock 
is a predictable classical behavior. In the range of (p where quantum fluctuations are 
important, (p is not suitable for this role, and the only remaining variable to be used as a 
clock is a. 

Another possibility is to abandon the principle of mediocrity and invoke the ideal 
observers that we used to define the distribution P(t), rather than physical observers, to 
calculate probabilities. In a recycling universe, the worldline of an ideal observer crosses 
an infinite number of inflating and thermalized regions. The probabilities for different 
types of thermalized regions can then be defined as relative frequencies at which these 
regions are encountered along the worldline. This definition is obviously gauge- invariant. 
In the model of Section V.B, it gives 



34 



PlM = \ Jl/J2 



(86) 



The gauge- invariance of (|^) is easily verified from Eqs. (^5|) , (^Tf) , (|50D , (|66|) . In tliis ap- 
proach, the most probable thermalized regions may turn out to be unsuitable for life, but 
this can be easily fixed by defining appropriate conditional probabilities. 

As mentioned at the end of Section III, recycling may not be complete in models 
where there is a "Planck boundary" in the diffusion regime or where some of the minima 
of the effective potential have vanishing or negative effective cosmological constant. In 
this case, the worldlines of all ideal observers (except a set of measure zero) have a finite 
length, and a natural extension of Eq. (^) is 



This defines the probabilities as being proportional to the total number of encounters 
for a given type of region, averaged over all observers. The result depends on the initial 
distribution P(0) at t = 0. Assuming that in this type of models the universe must have 
a beginning, and that it can be described by quantum cosmology, this initial distribution 
can be determined from the wave function of the universe. 

Although the definition of probabilities in this approach is gauge-invariant, it is not 
quite satisfactory. The ideal observers have very little to do with real physical observers, 
and it is hard to justify why the likelihood of various observations made by our civilization 
should be related to an ensemble of such ideal observers. 

We have to conclude that, despite some effort, none of the approaches suggested so 
far appears to be particularly compelling. It may turn out that, after all, an eternally 
inflating universe does not admit a uniquely defined probability distribution for the con- 
stants of nature. If so, this does not necessarily mean that all possible sets of constants 
consistent with our existence are equally likely. Although the ratio Pi/p2 may depend 
on the choice of cutoff procedure, it is conceivable that in some cases pi/p2 ^ 1 for all 



_ dt\J,{t)\ 

P2 J^dt\J2{t)\' 



(87) 
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reasonable choices. We would then "predict" that 1 is much more likely than 2. It is 
possible that we will have to restrict ourselves to such "stable" predictions, which are 
insensitive to the choice of cutoff . 



VII. CONCLUSIONS 

We have shown that the picture of the superlarge-scale structure of the inflationary 
universe is significantly modified by quantum fluctuations which bring parts of already 
thermalized regions back to the false vacuum, a process which we call recycling. 

In particular, the question of whether or not the Universe had a beginning is re- 
opened. Ignoring recycling, and under certain rather general conditions, it has been shown 
in Refs. that inflationary models require a beginning in time. In the recycling pic- 

ture, this conclusion does not apply because it is only necessary that any false vacuum 
region has a beginning in time. One can therefore imagine a nested structure where all 
false vacuum regions are just bubbles which nucleated inside preexisting true vacuum 
bubbles, which in turn nucleated inside false vacuum bubbles and so on. Whether or not 
this pattern can be continued to the infinite past is an interesting open question. 

We have extended the standard formalism of stochastic inflation to study 

the probability distribution of phases in which co-moving observers (or co-moving volume) 
find themselves in the recycling universe. Instead of the Fokker-Planck equation, one now 
has a more general master equation. We found that, in the case of complete recycling, all 
solutions of this equation rapidly approach a stationary form. This is in contrast to the 
standard case, where the probability disribution decreases exponentially with time, due 
to the loss of co-moving volume at thermalization. 

In the absence of recycling, the fractal dimension of the false vacuum co-moving 
volume at large times is lower than three. Including recycling, this dimension is d = 3 
because co-moving volume in false vacuum is continually replenished by nucleating false 
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vacuum bubbles. The universe ends up in a highly convoluted state, where the fractal 
dimension of any connected false or true vacuum region is lower than three. 

Finally, we have considered the question of "making predictions" for the constants 
of nature in the context of a recycling universe. The Principle of Mediocrity has been 
invoked in the past in order to obtain probability distributions for the constants. For the 
case of finite inflation the procedure is unambiguous: the probability is proportional to 
the number of civilizations that observed a given set of constants in the entire history of 
the universe. In the case of eternal inflation, the Principle is not so easy to implement, 
because the number of thermalized regions with given values of the constants is infinite 
(even in a finite co- moving region), and one has to introduce a regulator. If one simply 
counts all civilizations below some cut-off time, then the result depends strongly on the 
choice of time variable A gauge invariant cut-off prescription (i.e. one which does not 
depend on the time variable) was introduced in but this prescription is not unique 
||28|| . To make matters more complicated, the methods discussed so far cannot be directly 
applied to a recycling universe. We have considered some generalizations and alternative 
approaches, but none of them is particularly compelling. 

Therefore, it seems that while the Principle of Mediocrity may offer some valid guid- 
ance, it may not be sufficient to unambiguously determine the probability distribution for 
the constants of nature. If this turns out to be the final answer, then "predictions" would 
only be possible in those cases when all reasonable cut-off methods or implementations 
of the Principle yield a similar answer. 
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FIGURE CAPTIONS 

Fig. 1- True vacuum bubbles (white) nucleating in false vacuum (black). The shaded rings 
represent slow roll regions (external ring) and matter or radiation dominated regions 
(internal ring). 

Fig. 2- Self interaction potential for the "tunneling" scalar field. The energy densities in 
false and true vacua, denoted as p/ and pt, act as an effective cosmological constant. 

Fig. 3- The potential for the scalar field in the one field model of open inflation. The 
sharp barrier followed by a flat plateau gives this potential a somewhat unnatural 
appearance. 

Fig. 4- A false vacuum bubble (black) nucleating in true vacuum (white). Regions of slow 
roll and of matter and radiation domination surrounding the bubble are indicated. 

Fig. 5- The infiaton potential for the case of new inflation. 

Fig. 6- De Sitter space can be viewed as a hyperboloid embedded in a 5-dimensional 
Minkowski space. Here we represent a section of this hyperboloid, along the plane 
T, w in the embedding space. The region w + t > of the hyperboloid can be 
covered with flat FRW coordinates x, t. A section w + t = const, corresponds to 
t — const. 

Fig. 7- The Coleman-de Luccia instanton for bubble nucleation can be obtained by match- 
ing two 4-spheres of radii HJ^ and Hf^. The two spheres join at a 3-sphere which 
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represents the Euclideanized worldsheet of the domain wall. 

Fig. 8- The spacetime representing a false vacuum bubble in true vacuum. 

Fig. 9- Conformal diagram of the false vacuum bubble spacetime, for the Lee and Weinberg 
model. 

Fig. 10- Same as in Fig. 9 but for the case of a one field model of open inflation. 
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